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Shaouly Bar-Shalom 
Theoretical Physics Group, Rafael, Haifa 31021, Israel 

Gad Eilam and Ya-Dong Yang 
Physics Department, Technion- Israel Institute of Technology, Haifa 32000, Israel 

We study the pure penguin decays B — > 4>ir and B — » (fxf). Using QCD factorization, we find 
B(B ± -> (j)^) = 2.0t£? x W" 8 - For the pure penguin annihilation process B° — > cj><j>, analyzed 
here for the first time, B(B° — > (jxf>) = 2.Utg'g x 10 -9 . The smallness of these decays in the Standard 
Model makes them sensitive probes for new physics. From the upper limit of B — > 4>tt, we find 
constraints on R parity violating couplings, \\" i23^" i2i\ < 6 x 10~ 5 , ^23-^21! < 4 x 10~ 4 and 
|A-3 2 A-i 2 | < 4 x HT 4 for i = 1,2,3. Our new bounds on | A"i23 A"i2i | are one order of magnitude 
stronger than before. Within the available upper bounds for | A"i23 A"i2i | , I A^23 A^ 2 i I and | Ai 32 A^ 12 1 , 
we find that B(B — > cfxf>) could be enhanced to 10 -8 ~ 10~ 7 . Experimental searches for these decays 
are strongly urged. 
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Charmless two-body nonleptonic decays of B mesons 
provide tests for the Standard Model (SM) at both tree 
and loop levels. They also test hadronic physics and 
probe possible flavor physics beyond the SM. In past 
years, we have witnessed considerable progress in studies 
of these decays. Many such processes have been mea- 
sured or upper-limited Q. Theoretically, QCD factoriza- 
tion in which non-factorizable effects are calculable was 
presented Q , while there is also progress in perturbative 
QCD approaches B. 



In this letter, we will use QCD factorization to study 
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dominated by electroweak penguins, were studied by em- 
ploying naive factorization. By naive factorization one 
means that the hadronic matrix elements of the relevant 
four-quark operators are factorized into the product of 
hadronic matrix elements of two quark currents that are 
described by form factors and decay constants. In con- 
trast to color allowed processes, where naive factoriza- 
tion works reasonably well, this assumption is question- 
able for penguin processes. QCD factorization Q can 
be used to calculate the non-factorizable diagrams. We 
will use this framework to improve the theoretical pre- 
dictions for B '^ — > <frn 0,± . It is interesting to note that 
_B 0,± — > (fnr 0,± do not receive annihilation contribution, 
while B° — > (fxj) is a pure penguin annihilation process. 
To the best of our knowledge, there is no realistic theo- 
retical study of B° — > <fxj>. With respect to the topology 
of non-factorizable penguin diagrams for charmless B de- 
cays, the decay B° — > (fxj) is of interest. It can give us 
insight into the strength of the annihilation topology in 
nonleptonic charmless B decays which is still in dispute 
|3~|. Experimentally, B° — > (fxj) is relatively easy to 

x 10~ 9 and 
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|^|, Experimentally, B° — > (jxj) is re 
identify. We find B(B° -> (fxf)) = 2.lj 
B(B ± — ► (j)ir^~) = 2.0±g f x 1CT 8 in the SM. The small- 
ness of the SM predictions for these decays makes them 
sensitive probes for flavor physics beyond the SM. We 
use the recent Babar upper limits B(B°' ± — > </>7r 0,± ) < 



1.6 x 1CT 6 , at 90% CL [§, to obtain limits on the rele- 
vant R Parity Violating (RPV) couplings. We then use 
these limits to deduce the maximal possible enhancement 
of B(B° — ► (fxf)) in RPV supersymmetry. Note that cur- 
rently B(B° -> (jxj)) < 1.2 x 1CT 5 at 90% CL §. 

The experimental upper limits could be improved in 
BaBar and Belle. Measurement of any of these decays 
with B )i 10~ 7 will serve as an evidence for new physics. 

In the SM, the relevant QCD corrected Hamiltonian is 
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The operators in 7i e // relevant for b — > dss are given 
in ||, where at the scale \i — nib, C3 = 0.014, C4 = 
-0.035, C 5 = 0.009, C 6 = -0.041, C 7 = -0.002/137, 
C 8 = 0.054/137, C 9 = -1.292/137, C 10 = 0.262/137. 
Using the effective Hamiltonian and naive factorization 
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and A(B° -> 0tt°) = -^A(B~ -> <jyK~) with a, = 
Ci + Ci+\/N c . The contributions of strong penguin oper- 
ators arising from the evolution from fi = Mw to fi = mb 
is very small due to the cancellations between them: 
C3(m b ) ~ -C 4 (m b )/3 and C 5 (m fc ) ~ -C 6 (m b )/3. Ob- 
viously the amplitude is dominated by electroweak pen- 
guin. Using / = 254 MeV, \V td \ = 0.008 7V C = 3, and 
the form factor Ff^O) = 0.28 §, ||, y, O, we get 
B(B ± -> tjm*) = 2B(B° -> 0tt°) = 2.9 x 10"^ 

In the above calculations, non-factorizable contribu- 
tions are neglected. However, this neglect is question- 
able for penguin dominated B — > (/jtt. The leading 
non-factorizable diagrams in Fig.l should be taken into 
account. To this end, we employ the QCD factoriza- 
tion framework 0], which incorporates important the- 
oretical aspects of QCD like color transparency, heavy 
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quark limit and hard-scattering, and allows us to calcu- 
late non-factorizable contributions systematically. In this 
framework, non-factorizable contributions to B~ — > n~ '<j> 
can be obtained by calculating the diagrams in Fig.l. 
To leading twist and leading power, the amplitude for 
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FIG. 1: Non-factorizable diagrams for B — * 4>tv 
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and A'(B° -► 0tt o ) = ^A'(£f~ -> 07r _ ). The a s term 
is the non-factorizable contribution with 

F«i = -12 In— - 18 + V + S, 



(4) 

v = J q du$4,(u) (s 1 ^ 2 " Inn - 3ztt^ , (5) 
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where £ = 1 + /Mb is the momentum fraction carried by 
the spectator quark in the J3 meson. The $'s are the lead- 
ing twist light-cone distribution amplitudes of tt, </> and 
B mesons. They describe the long-distance QCD dynam- 
ics of the matrix elements of quarks and mesons, which is 



factorized out from the perturbative short-distance inter- 
actions in the hard scattering kernels. These distribution 
amplitudes can be found in p3j |l4l |±5|]. In our calcula- 



tion, we use the model proposed in 
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where A is the momentum scale of the light degrees of 
freedom in the B and taken to 350 MeV. To show model 
dependence of our prediction, we vary A from 150 MeV 



to 550 MeV. we get B(B d 



--) = 2B(B° </m°) = 



2.0^;? x 10~ 8 . From Eq.| we see that non-factorization 
is dominated by strong penguin due to the absence of 
Cg. We also note that non-factorizable contributions 
dominate these decays and there is no isospin symmetry 
breaking because annihilation contributions are absent. 

The decay B° — ► <fxp is also of interest to study. Firstly 
it is a pure penguin process. Secondly it is a pure anni- 
hilation and thirdly its experimental signature is very 
clean. By naive factorization, the amplitude for this de- 
cay mode is 
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This amplitude vanishes for m s — > 0. The a s order ma- 
trix (4>(f>\s $b (1— tb)s|0) also vanishes due to the cancel- 
lation between the amplitudes of Fig. 2. (c) and (d). Non- 
factorizable contributions can be obtained by calculating 
the amplitudes of Fig. 2. (a) and Fig.2.(b). They are 
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where we set \ih to be the average virtuality of the time- like gluon, fj,h = Mb/2. In Eq|T(], fc 2 and D b ^ d are the 
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FIG. 2: (a) and (b) are non-factorizable diagrams for B° — > (f><j> decays, (c) and (d) are factorizable diagrams at a„ order. 



virtualities of gluon, b and d quark propagators, respec- 
tively. As in |L6| ], we meet end point divergence when 
U 



0. Instead of a cut-off treatment Jl6| , 



we use an 



effective gluon propagator |Tj 
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Typically m g = 500 ± 200 MeV, A = A QCD = 300 MeV. 
Our use of this gluon propagator instead of imposing a 
cut-off, is supported by lattice |jTs| , and field theoretical 
solutions ]l9| ] which indicate that the gluon propagator 
is not divergent as fast as *r. Finally we get B(B° — > 
H>) = 2.1+i;l x 10- 9 . 

Potentially, this decay may be enhanced by rescatter- 
ing B — > rf'^rf-'} — > 4>4> or by uj — <f> through the channel 
B° — > uj(f> — > (fxj). However, n' and r\ contributions are 
almost completely canceled [f20|, and <fi is nearly a pure ss 
state, so the mixing mechanism is also negligible. Fur- 
thermore in the language of QCD factorization frame- 
work, such kinds of soft final state interactions are sub- 
leading and suppressed by power of O ( Aqc D/mb) , al- 
though it is hard to be calculated reliably. Lastly, strong 
interaction annihilation is negligible since at least two 
gluons should be exchanged. Thus, any unexpected large 
branching ratio observed, will indicate new physics. 

As an example for new physics, we will discuss the 
effects of the trilinear A' and A" terms in the RPV su- 
perpotential p2, 23, 23 on the process b — + dss. We 



are therefore interested in 
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where a, b are 577(2) indices, i, j, k are generation indices, 
a, /?, 7 are SU(3) color indices and c denotes charge con- 
jugation. The L (Q) are the lepton (quark) SU{2) dou- 
blet superfields, and U (D) are the up- (down-) quark 
SU{2) singlet superfields. Then we have 
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From , we get the effective Hamiltonian for b — > dss 
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The coefficients 



r\ 4 /P° and rj 8 /^° are due to running from the sfermion 
mass scale mi (100 GeV assumed) down to the nib scale. 

We can now write down the contributions of Ti.^ to 
B~ — ► </>7r~ and B° — > 00 decays 
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In the numerical results, we assume that only one 
sfermion contributes at a time and that they all have 
a mass of 100 GeV. The uncertainties of the theoreti- 
cal predictions, due mainly to the B meson distribution 
function, are displayed as thickness of curves in Fig. 3 

Our results for the RPV contributions to B — > <\>-k are 
summarized in Fig. 3. From the BaBar upper limit || 
B(B°' ± -> cw ^) < 1.6 x 10" 6 , we obtain the following 
constraints (90%CL) 

lA^A^I < 6xHT 5 (^) 2 , (17) 
142^1 < 4xl0- 4 (^) 2 , (18) 

ia; 21 a^ 3 i < 4xio ~ 4 (i§f) 2 • ( 19 ) 

We note that our constraints on A i23 A i2 * are more than 
one order of magnitude stronger than the limits obtained 
recently [^5| . For A i32 A i * 2 and A i21 A i23 , our bounds 
are comparable with the present upper limits |25|, p6| . 
Within the available upper bounds for these couplings, 
X'- 23 X'-2i RPV terms could enhance B(B° -> <t>4>) to 10~ 8 , 
while A- 23 A- 21 and A- 32 A-^ 2 RPV terms could enhance 



B(B° — > 4>(j>) to 10 7 which may be measurable at Belle 
and Babar. 

In summary, we have studied the pure penguin pro- 
cesses _B ±! ° — > (fnr^' and B a — > (fxp by using QCD 
factorization for the hadronic dynamics. We estimate 
that in the SM B(B~ -► <jm~) = 2.0±g f x 10~ 8 and 
B(B° -> <fxj>) = 2.1±J;1 x 10~ 9 . The smallness of these 
decays in the SM makes them sensitive probes of fla- 
vor physics beyond the SM. Using the BaBar result 
B(B~ — > 07r~) < 1.6 x 10~ 6 , we have obtained new 
bounds on some products of RPV coupling constants. In 
the case of A" 23 A'/-i* 2 , our limits are better than previous 
bounds. Given the available bounds on A i23 A i2 \, A^ 23 A£ 21 
and A- 32 A-^ 2 , the decay B° — ► <jxf) could be enhanced 
to 10~ 8 ~ 10~ 7 . Due to the clear signatures of 4> and 
n^ 1 , the experimental sensitivity of these decay modes is 
high. Babar and Belle could reach very low upper lim- 
its on these decays if not measured. Searches for these 
decays are strongly urged. 
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FIG. 3: The branching ratio of B~ — * (f)n~ as a function 
of the RPV couplings |A" 23 A"i* 2 | (upper curve), |A- 23 A- 21 | and 
I Ai 3 2 Aii2 1 (lower curve) respectively. The thickness of curves 
represent our theoretical uncertainties. The horizontal lines 
are the upper limits and the SM prediction as labeled respec- 
tively. The thicknesses of the curves and the line labelled as 
SM are theoretical uncertainties. 



